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We study a one-dimensional lattice model of interacting spinless fermions. This model 
is integrable for both periodic and open boundary conditions, the latter case includes the 
presence of Grassmann valued non-diagonal boundary fields breaking the bulk U(l) sym- 
metry of the model. Starting from the embedding of this model into a graded Yang-Baxter 
algebra an infinite hierarchy of comuting transfer matrices is constructed by means of a 

o- 

^vq ■ fusion procedure. For certain values of the coupling constant related to anisotropics of the 

" underlying vertex model taken at roots of unity this hierarchy is shown to truncate giving 

o: 

a finite set of functional equations for the spectrum of the transfer matrices. For generic 



(N 



> 
oo 

(N 



coupling constants the spectral problem is formulated in terms of a TQ-equation which can 
be solved by Bethe ansatz methods for periodic and diagonal open boundary conditions. 
Possible approaches for the solution of the model with generic non-diagonal boundary fields 



are discussed. 

•5: 



INTRODUCTION 



The small polaron modelprovides an effective description of the behaviour of an additional 
electron in a polar crystal [ll, [3] - In one spatial dimension this lattice system of interacting spinless 



:*ermions can be constructed within the framework of the Quantum Inverse Scattering Method 

n nn 

3] allowing to compute the excitation spectrum by Bethe ansatz techniques, see e.g. [4|, |5j|. By 
means of a graded generalization of Sklyanin's reflection algebra [9] it was possible to provide 
the small polaron model with open boundary conditions while keeping its integrability intact. 
These integrable boundary conditions are encoded in c-number valued 2 x 2-matrix solutions to 
the reflection equations [1CH12 ] . 



Diagonal boundary matrices correspond to boundary chemical potentials in the Hamiltonian. 
In this case the spectrum of the open small polaron model can be studied again using Bethe ansatz 



methods 



13l4l5||. As a consequence of the underlying grading the general non-diagonal solutions 



to the reflection equations are super matrices 



Iff]. The resulting additional boundary terms do 



not conserve particle number and have anti-commuting scalars, i.e. odd Grassmann numbers, as 
amplitudes. The fact that the U(l) symmetry of the model is broken implies that there is no 
simple eigenstate (e.g. the Fock vacuum) of the model which can be used as a reference state for 



2 



the algebraic Bethe ansatz. This situation is, in fact, very similar to the case of non-diagonal 
boundary magnetic fields in the closely related spin-1/2 Heisenberg chains: in all approaches used 
so far the solution of the spectral problem is limited to small finite systems, with constraints 



between the boundary fields at the two ends of t 
their usefulness in the thermodynamic limit 17 



re chain or restrictions on the anisotropy limiting 



25]. 



In a previous publication [26[ we have investigated the applicability of Bethe ansatz methods 
in the simpler case of a model of free fermions with similar open boundary conditions. We found 
that for a certain class of non-diagonal boundary super matrices, a unitary transformation on 
the auxiliary space allowed for an exact solution of the free fermion model. Furthermore, the 
functional equations obtained there could be easily generalized to describe the spectrum of the 
model for arbitrary non-diagonal boundary fields. Unfortunatley, this approach can not be applied 
to the small polaron model. 

Following ideas 



[13,120] 

of the spin-1/2 XXZ chain before we adapt the fusion 
procedure for the transfer matrix of the quantum chain to the graded case of the small polaron 
model. Assuming the existence of a certain limit we obtain a system of functional TQ-equations 
for the eigenvalues of the transfer matrix which are shown to coincide with the result obtained from 
the algebraic Bethe ansatz for both periodic and diagonal open boundary conditions. For special 
values of the interaction parameter related to roots of unity of the anisotropy parameter we derive 
truncation identities for the fusion of the relevant objects, in particular the transfer matrices. Using 
these identities the TQ-equation reduces to a a set of relations between finitely many quantities. 



II. THE SMALL POLARON AS A FUNDAMENTAL INTEGRABLE MODEL 

Some materials exhibit a strong electron-phonon coupling that considerably reduces the mobility 
of electrons within the conduction band. This interaction may be regarded as an increase of the 
electron's effective mass, thus giving raise to quasi-particles called polarons. If the electron is 
essentially trapped at a single lattice-site the corresponding quasi-particle is said to be a small 
polaron. In this case, electron transport occurs either by thermally activated hopping (at high 
temperatures) or by tunneling (at low temperatures). 

In the case of periodic boundary conditions (PBC) the iV-site small polaron model is charac- 
terized by the Hamiltonian 

N 

jjPBc _ w ith Hn,n+i = Hn,i (2.1) 

3=1 
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with a Hamiltonian density Hjj + i defined as 

where c£ and c fc label the creation resp. annihilation operators of spinless fermions at site k, which 
are subject to the anticommutation relations [c£,c fe ] + = 5%. Moreover, it is convenient to define 
number operators n k = cic. = 1 — n, ■ In this context, the parameters t and V may be interpreted 
as hopping amplitude and density-density interaction strength respectively. 



A. Construction within the QISM framework 

The small polaron model can be associated to a graded six-vertex model with anisotropy rj and 
i?-matrix 



R(u) = 



sin(2?7) 



sin(u + 2rj) 







\ 



V 



sin(u) sin(2?7) 

sin(2?7) sin(u) 
-sm(u + 2ri)J 



R(u) a solution to the Yang-Baxter Equation (YBE) 



(2.3) 



R± 2 (u - v)Ri 3 (u)R 23 (v) = R 23 (v)Ri 3 (u)R 12 (u - v) 



(2.4) 



and enjoys several useful properties, such as 



P-symmetry 



R 2 i(u) = V\ 2 R\ 2 {u)V\ 2 = R\ 2 {u) 



(2.5a) 



T-symmetry 



Rf 2 ^(u) = R^(u) = R 21 (u) 



(2.5b) 



regularity 



#i 2 (0) = V 12 



(2.5c) 



• unitarity 



Ri 2 {u)R 21 (-u) = C(«) 



(2.5d) 
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where the scalar function £(u) is given by 

C(u) = g(u)g(-u) and g(u) = . . 

sin (277 J 

Unitarity of an i?-matrix is of course a direct consequence of its regularity. 
crossing symmetry 

R% (-u - 477X 1 (u) = ((u + 2r?) (2.5e) 

periodicity 

R 12 (u + vr) = -cr| i? 12 (u) of = -of i?i 2 («) of (2.5f) 
The periodicity R(u + 27r) = R(u) is obvious from definition 



The operations of partial super transposition st a and inverse partial super transposition ist a as 
well as the graded permutation operator V a b and the notion of super tensor product structures are 
explained in appendix [A] Unless stated otherwise, all embeddings are to be understood in a graded 
sense, that is into a super tensor product structure. Considering the Yang-Baxter Algebra (YBA) 

R l2 (u - v) T x {u) T 2 {v) = T 2 (v) T x {u) R 12 (u - v) (2.6) 

this means that Ti(u) = T(u) ® s 1 and T 2 (v) = 1 (g> s T(v). 

The small polaron model constructed here is fundamental, i.e. the Lax-operators Lj(u), being 
local solutions to (|2.6p are just graded embeddings of the above i?-Matrix (|2.3p . 



1 / sm(u)ri: + sin(n + 2rf)n • sin(2r/)ct \ 

Lj(u) = . I " 

sm{zr]) \ sm(2r])Cj sin(it)ra- — sin(« + 2rj)n^ J 

As a consequence of the YBA's co-multiplication property, a specific global representation, the 
so-called monodromy matrix, can be constructed as a product of Lax-operators taken in auxiliary 
space, 

T{u) = L N {u) ■ . . . ■ L 2 (u) ■ L^u) (2.8) 
and gives rise to a family of commuting (super) transfer matrices 

t(u) = str { T(u) } [t(u),t(v)] =0 Vu,d£C, (2.9) 

where str { • } denotes the supertrace defined in appendix [A] In particular, the PBC hamiltonian 



2.ip with t = 1 and V = — cos(2?7) is among these commuting operators, 

(2.10) 



H PBC = - sinter?)— In r(n) 
du 



u=0 
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B. Asymptotic behaviour of the PBC transfer matrix 



By construction the monodromy matrix (and similarly the transfer matrix) is a Laurent poly- 
nomial in z = e m , i.e. T(u) = J2k=-N TkZ k ■ For z — > oo the Lax-operators ()2.7p are 



LAu) 



n- + e 2 " 7 n- 







n- — e 2lv n- 



2isin(2?7) y (! 

and consequently the asymptotic behaviour of the (super) transfer matrix is given by 



(2.11) 



t(u) 



2i sin(2ry) 



N 



jNn 



N N 

n(e- i s+e i s)-n( e " i s- ei s 

3=1 3=1 



(2.12) 



As the leading term comprises only diagonal operators the first order contributions to the transfer 
matrix eigenvalues Am(u) can easily be determined and are found to depend on the total number 
of particles M, 

N 



A M (u) « e 



iuN 



e i7Vj7 e -iM2r; 



M -iNrj JM2ri 



■ire 



gi2?7 _ g— i2r) 

This result will be used to fix the degree of the Q-functions in section IIVI 



) • 



(2.13) 



III. FUSION OF THE i?-MATRIX IN AUXILIARY SPACE 

Given an R-matrix as solution to the YBE (|2.4p . the fusion procedure allows for the construction 
of larger R-matrices as solutions to corresponding Yang-Baxter equations, where larger refers to the 
dimensionality of the auxiliary space involved. All that fusion requires is a pair of complementary 
orthogonal 1 projectors P-jt an d P\2 such that for a specific value of p G C the following triangularity 
condition holds for arbitrary spectral parameters u G C, 

P u R 13 (u) R 23 (u + p) P+=0 . (3.1) 

By virtue of this condition, it can be shown that the fused R-Matrix, defined by 

fl(12) 3 («) = Pl 2 Rl3(u)R23(u + P )P+ . (3.2) 

satisfies the corresponding Yang-Baxter equation 

R(12)3(u ~ V) R(i2) 4 (u) R3i(v) = R 3 a(v) #(12)4 0) #(12)30 ~ ") ■ (3-3) 
1 As usual, orthogonal means P^~ 2 P^ 2 — whereas complementary refers to the property P^ 2 + P^ 2 — 1. 
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It is easily found that the small polaron R- matrix (|2.3p has two distinct singularities at u = ±2ry, 



det{R(u)} 



sin(n — 2rj) ( sin(ii + 2rj) 



0. 



(3.4) 



(3.5) 



sin(2?7) y sin(2?7) 

At u = —2rj the R-Matrix gives rise to a projector onto a one-dimensional subspace, 

/o (A 
1-10 
0-110 
\0 0/ 

However, unlike in the case of the Heisenberg spin chain, the orthogonal projector P + onto the 
complementary three-dimensional subspace cannot be obtained from the R-matrix at the second 
singularity, 



p + ^t-p-^ ~R(2n) . 



(3.6) 



Using this projector, fusion of two small polaron R-matrices in the auxiliary space can be achieved 
by means of (|3.2|) with p = 2r], 



R(i2)3(u) = P+Rn(u)R 23 (u + 2rj)P+ 



(3.7) 



The resulting object R(i2)3( u ) 1S an 8x8-matrix of rank 6 and may therefore be effectively reduced 
to a 6x6-matrix -R<i2>3(^) acting on a three-dimensional auxiliary space V^i2> and on a two- 
dimensional quantum space V3. Changing from the BFFB-gr&ded 2 canonical basis 



B = {ex, e 2 , e 3 , e A } BFFB = {|0) ® |0), |0) ® |1), |1) ® |0), |1) ® |1)} BFFB 
to the projectors' BFBF-graded singlet/triplet-eigenbasis 

B± = {/l> /2; /3; f^BFBF = 7= 3 : 64, .- 3 }bfbf • 



the matrix -R(i2)3(^) gains the advantageous shape 



^<12>3(w) 


\ 









= (/if [i2(i 2)3 («)] /j 


V 


0j 





(3.8) 



(3.9) 



(3.10) 



This notation is explained in appendix 1X1 
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where -R<i2>3(' u ) is the only non- vanishing block. Explicitly one finds, 



#<12>3W oc 



/ 

' 2 sin(tt+4?)) 











1 





2 sin(u) 


v / 2sin(4»)) 











2-\/2sm(2r?) 


2sin(«+2r?) 














-2sin(u+2r?) 


^v^sin^?;) 











v^sin^) 


2 sin(u) 





V o 











2 sin(it+4?i) J 



(3.11) 



A. General construction of higher fused i?-matrices 

In general, higher fused R-matrices can be constructed employing the projection operators 

(3.12) 



<r6<S n 



Here a runs through all the elements of the permutation group S n and P a is the permutation 
operator corresponding to a. Now the higher fused R-matrices are obtained as 



R(i...n) q {u) = P+. n Ri q (u) R 2q (u + 2 V ) ... R nq {u + [n - 1] • 2r ? )P+ „ 



(3.13) 



Just as for the first fusion step, it is convenient to apply a similarity transformation A(i ---n ) into 
the eigenbasis 3 of the projection operators, 



A( 1 ...n)R(l... n ) q (u)A (1 1 n) = 





f 


\ 


n) ~ 









V 





(3.14) 



The first few (n = 1,2,3,4) transformation matrices -A(x... n ) are explicitly given in appendix [Ej By 
construction, all matrix elements of (|3,14p . except for those in the upper left 2(n + 1) x 2(n + 1) 
block, vanish. This block is referred to as the fused R-matrix i?<i... n >g("u). As shown in table HI 
its fused auxiliary space has alternating gradation (bosonic, fermionic, . . . ). 

auxiliary space: <12> <123> <1234> <12345> 
grading: BFB BFBF BFBFB BFBFBF 



3 Since the projectors here are just the same as for the XXZ Heisenberg spin chain, the respective transformation is 
simply given by the matrix of Clebsch-Gordan coefficients. 
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TABLE I: Gradation of the fused auxiliary spaces in the projector eigenbasis. 

The periodicity property (|2.5f|) carries over to the fused ^-matrices, 
#<l...n»g(M + 7r) = (-l) n CJ< <n>> l?«i...„» g ()i) a\ 
with o"^- n> being defined through 



(n) 



n 

fc=i 



o-| and ^4(12... n) of„) ^(iL.n) 



\ 



(3.15) 



(3.16) 



B. Fusion hierarchy for super transfer matrices 

Since, by construction, the fused i2-matrices again satisfy the YBE they can be used to establish 
further families of commuting operators as supertraces of fused monodromy matrices, 



(I2...n) 



n){u) 



P 



12. ..n 



R(12...n)q N (u) R(l2...n)q 2 ( u ) R(12...n)qi( u ) P l 



12...n 



Pi + 2...n P(i2...n-i)(«) T„(u + [n - 1] • 2r?) P+... n . 



4 



(12. ..n) i (12...n) 



:n) a; 





/ 


\ 




?<i2...n>('") 




.n) 









V 





(3.17) 



(3.18) 



Indeed it is found that the (super) transfer matrices obtained from any fusion level n, 

r (n) (u) = str (12 n+1) { T {12 ... n+ i)(n) } = str <<12 n+1>> { r <12 ...n+i>(«) } , (3.19) 

commute with the transfer matrices of any other fusion level m, i.e. [r (n) (n), r (m) (t>)] = for all 
u, v E C and arbitrary n, m £ No- A most interesting fact is, that these fused transfer matrices 
obey certain functional relations, known as fusion hierarchy. For the periodic boundary case, the 
fusion hierarchy reads, 

r (0 \u) r (0) (n + 2r/) = r (1) (n) + 5{u) 

t w {u) r (0) (n + 4t/) = r (2) (n) + 5{u + 2r/)r (0) (n) 

r (2) (u) r (0) (u + 6r/) = r (3) (u) + <5(u + 4r/)r (1) («) (3.20) 



r (n) («) r (0) ( u + [ n + 1] • 2t/) = r (n+1) («) + *(« + n ■ 2r])T (n ' r) (u) , 
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where 5(u) = 5 {T(u)} labels the PBC super quantum determinant (SQD) defined in appendix [Dl 
In contrast to ungraded models, such as the XXZ Heisenberg spin chain, this quantum determinant 
is not proportional to the identity. 



IV. TQ-EQUATIONS FOR PBC 

After applying a shift u — > u — [n + 1] • 2rj the PBC fusion hierarchy (|3.20p reads 

rW(«- [n + 1] -2rf) t (0) {u) = r (n+1) (u - [n + 1] • 2rj) + 5(u - 2rf) r^" 1 )^ - [n + 1] • 2 V ) . (4.1) 

As all operators in this equation mutually commute, it may equally well be read as an equation 
for the eigenvalues A^ n \u) of the fused super transfer matrices. With A(u) = A^°\u) this yields 

A{U) ~ A(n)(u-[n + l]-2 V ) { ' ^ 1 ' AW (n- [n + 1] -27?) ( ' 

where M is the number of particles in the system, such that the sign (— 1) M depends on the parity 
of the corresponding eigenstate (bosonic/fermionic). This pecularity stems from the fact, that the 
PBC SQD (|D21aP can not simply be treated as a scalar function but rather as an operator that 
intersperses sign factors into the respective sectors. This may be illustrated by considering the 
fusion hierarchy (|4.1[) in a diagonal basis for chain length N = 1, 





B 
F 



(4.3) 



Introducing the functions 

Q W (u) = A< n) {u - [n + 1] • 2??) (4.4) 
the eigenvalues can be rewritten as 

Now factorize according to 

Q^=Xm{u)-T^{u).Q^{u) (4.6) 



where 



(4.7) 
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Provided the existence of the limit Q(u) = lim n _ 5 . OC) Q(")(u), this yields 

= / sin(n + 2r]) \ N Q(u -2rj) M ( sin(u) \ N Q(u + 2rj) 

{U) V sin(27?) J Q(u) [ 1 \sm(2ri)J Q(u) 

Due to the structure of the entries in the Lax-operators, the Q-functions factorize 

G 

Q(u) = Y\ sin(n — A^) 



(4.8) 



(4.9) 



where the integer G can be determined by considering the asymptotic behaviour of A(u). In the 
limit z = e lu — y oo the leading contribution to (|4.8p is 



A(«) 



ANu 



TV 



jNr, „-iG2 v 



_l)M e -iN v e iG2rn 



(4.10) 



gi 2r/ g— i 2rj 

such that consistency with (|2.13p immediately fixes G = M. The requirement for the eigenvalues 
A(u) to be analytic ultimately yields 



Res A ,(A) = * (f^+plYl Y[M*i-*i + *v) 



which are precisely the Bethe equations for this model 



sin(Aj) / sin(A^ — Xj + 2r\) 



(4.11) 



V. TRUNCATION OF THE PBC FUSION HIERARCHY 



In the case of the XXZ-model it has been observed, that for certain values of the anisotropy r] 
the fusion hierarchy repeats itself after a finite number of steps. The small polaron model shares 
this feature at values t] = T] p where 

tt/2 



Vp 



p + l 



(5.1) 



A. i?-matrix truncation 



The truncation identities for the R-matrices are found to be 



where 



Vp) 



(- 



\ 



■M p {u) a\ 



M p (u) {alY) 



TZf\u,7]) = fl<l...(p+l)»g(«) S <l...(p+1)> 

1/2 \ p sin([p + 1] u) 



M p (u) 



v sin(2 r] p ) J sin(2 rj p ) 
with the transformation matrices -B<i...n> explicitly given in appendix[E] up to n = 4. 



(5.2) 



(5.3) 
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B. Super transfer matrix truncation 

For periodic boundary conditions the S-transformed fused monodromy matrix T^ p \u,r]) of an 
TV-site model with quantum space H = V qi <8> s V q2 ® s . . . ® a V qN is defined as 

T^(u, V ) = n$>(u,r,) Kfl^v) ...K%>(u,t,) (5.4) 

= - B <l...(p+l)> - R <l...(p+l)>gjv( n ) ••• - R <l...(p+l)>gi( n ) - B <i...(p+i)> 
= £<l...(p+l)> T <l...(p+l)>( n ) jB <i...( p +i)> 
and due to the cyclic invariance of the supertrace it yields the exact same transfer matrix 

T {p) (u,ri) = str <<x {p+1)>> { r« 1 ... (p+1)> («) } = st r<<1 (p+1)>> { T (p) (u,ri) } . (5.5) 
At rj = rj p the truncation identity (|5.2p for R-matrices gives 
T^(u, Vp ) = 



(P-2), 



such that the truncation identity for the transfer matrices is found to be 

r^(u,r ]p ) = [-M p (u)] N (f[aA - (-l)P[M p (n)] N (UKn 



\i=l / \i=l 



(5.6) 



(5.7) 



VI. THE SMALL POLARON WITH OPEN BOUNDARY CONDITIONS 
A. Reflection algebras and boundary matrices 

The construction of integrable systems with open boundary conditions is based on representa- 
tions of the graded reflection algebra 

Ru(u - v)T{ (u)R 2 i(u + v)T 2 {v) 

(6.1) 

= T 2 {v)R\2{u + v)Ti {u)R 2 i{u - v) 
and the corresponding dual graded reflection algebra 

Rn(v - u)T^{uf^R 21 {-u -v- 4r/)T 2 + (v) ist2 

(6.2) 

= T 2 + (v) ist2 R 12 (-u -v- 4r ? )r i + (n) stl J R 21 (v - u). 
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The relation between R^u) and the conjugated R-matrix R a b(u) is explained in appendix [Bj c- 
number valued boundary matrices, compatible with the respective reflection equation are found to 
be 



K-(u) = uj- 
K + (u) = a; 4 



(6.3) 



sin(ii + a_ sin(2ti) 
{3_ sin(2u) — sin(n — if)_) 

sin(it + 2r? + ip + ) a + sin(2[ii + 2r?]) 
j3 + sin(2[u + 2rj\) sin(u + 2r? — ?/>_, 

with normalizations = w ± (??) defined by 

uj~(ri) = — — - and oj + (ri) = ; — — — r . (6.4) 

y " sin(^) 2003(27?) sin(^+) 

These matrices share the periodicity property of the i?-matrix, i.e. 

K*(u + tt) = -a z K*{u) o~ z . (6.5) 

Here the normalizations were chosen such that 

K-{0) = 1 and str { K+(0) } = 1 , (6.6) 

but apart from this, the two solutions are related via 

1 



K + (u) 



K-{-u-2q) a 1 



(e -> e) 



(6.7) 



2008(27?) 

where (0 — > ©) marks the replacements (a_, /3_, V-) ~~ ^ ( — a +, P+, — if>+)- In principle the 
parameters ip± are arbitrary even Grassmann numbers but their invertability requires them to 
have a non- vanishing complex part 4 . The remaining parameters a± and (3 ± are odd Grassmann 
numbers, being subject to the condition a± ■ /3± = 0. 

Given the monodromy matrix T(u) = Ln(u)Ln-i{u) . . . L\(u) it is possible to construct a 
further representation of the reflection algebra (|6.ip as 



, A(u) B[u) , 

T-{u)=T{u)K-{u)T{u)= ( | , (6.8) 

C(u) V(u) t 



4 Such an additive part, that contains no nilpotent generators, is sometimes called the body of a Grassmann number. 
It is to be distinguished from the soul of a Grassmann number, which contains only sums of products of nilpotent 
generators. 
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with T(u) being a shorthand notation for T ( 



-u 



f(u) = R^(-u) Ra h 



H)2 



-u 



Rqn(- u ) 



H5i 



< N (u) 



C(«) 



R 10 (u) R 2 q(u) . . . Rno(u) 



(6.9) 



N 



R i(u) R 2(u) ... Rqn{u) , 



resulting in an OBC super transfer matrix 

T(u) = sti {K + (u)T-(u)} 

Expanding t(u) around u = one obtains a Hamiltonian featuring the same bulk part (|2.2p 
corrsponding PBC Hamiltonian. Defining the shorthands M± = \ csc(2?7) csc(V>+) sin(2ry±-0. 
resulting OBC Hamiltonian 

iV-l 

jjOBC _ 



(6.10) 

as the 
A the 



+ 



A/"+ n N -AC % 



(6.11) 



+ csc(V>~) a- Cj — /?_ c{ + csc(t/> + ) a + — f3 + cjy 
is derived from the set of open boundary transfer matrices by 

d u t(u)\ u=0 = 2 H OBC + const. . 



(6.12) 



In the case of diagonal boundaries, i.e. a± = (5± = 0, Bethe equations can be derived using the 
algebraic Bethe ansatz. This allows for the computation of the transfer matrix eigenvalues and 
eigenvectors (cf. appendix Oresp. [3]). 



B. Asymptotic behaviour of the OBC transfer matrix 



Just as for the PBC case the highest order contributions to the (super) transfer matrix in the 
limit z = e m — > oo will be considered. In addition to the asymptotic behaviour ()2.1ip of the Lax 
operators Lj(u) it is found 



' M- o 



4sin(2r/) / n j + ^"' n } 
2iz 



+ o 



1 



K-(u) 
K + {u) 







( a\ 




~2T 






|-( 







o 



ur e 



2ir/ 



2i 



z 2 e 2ir > 



( a+\ 






H 








+ o 



(6.13a) 
(6.13b) 
(6.13c) 
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which yield the asymptotics of the OBC transfer matrix (|6.10p and of their eigenvalues Am(u) 
respectively 

N 



-1 



e 4ir ? - a + /3_) z 4 - e 2 ^.)^. + e 21 ^) + 0(; 



(6.14) 



A±(«)=±(-l) 



v 



e 4i " (/3+q_ - q+/3_) z 4 e i7V2r ? + 0(z 2 ) 



The eigenvalues A±(u) have been classified according to a parity which is determined by the 
(diagonal) operator controlling the asymptotics of t(u). 

Note that in the diagonal case, i.e. a± = J3± = 0, the 0{z) terms of the K matrices become 
relevant such that 

N 



t{u) = -(-1) 



jy U U 2 \r, J2 

4 6 



3 i(V'++V'-) 



V 



(6.15) 



Am(«) = "(-1)^ e 2i ^ 2 ^W + +^-) e 4i(JV-M)^ + e -i(V + +^-) e 4iM^ + 

Here, as in the case of periodic boundary conditions, the asymptotic behaviour of the transfer 
matrix eigenvalues can be related to the (conserved) total number M of particles in the state. 



C. Fusion of the boundary matrices 

For the sake of readability it is convenient to define the following ordered product of i?-matrices, 

i 

Rf^u) = Yl R k , i+l (2u + [i + k - 1] • 2r?) , (6.16) 
fc=i 

such that the fused K~ boundary matrices may be written as, 

Tn-l 



K (12...n)( U ) = P l2...n 



Y[K7( u +[i-l]-2r,) R: 



string 



U 



i=l 



K-( U + [n-l}-2 V ) P+.. 



.4 



(12.. .n) K { l2...n)( U ) A (12...n) = 



where K, 



(12. ..n) 



U 



( 


\ 


K «l...n^( U ) 









\ 


'•■ / 


.. n> (u) being 


the only 



(6.17) 



dimensions (n + 1) x (n + 1). There is a useful relation between the fused K - and i^ + -matrices 



15 



that stems from (|6.7p . 

K (12...n)^ 



1 



2cos(2r/) 



K 



(n...21)(- U ~ 71 ■ 2r l) °\n) 



( 



\l2...n) K^ n) (u) A { ^ n) 



K 



<l...n> 



(u) 



\ 



(6.18) 



/ 



and defines the (n + 1) x (n + l)-matrix i^J-i n >(w) in the obvious way, where <r?^ was defined in 
(|3.16p . Note that the order of all spaces in K7 21 % is inverted. Thus, by changing the space labels 
according to i — > n + 1 — i the fused right boundary matrix may explicitly be written as 

fn-l 

HK: +1 _ t iu+[n-i}.2 V )^(u) 



L (12...n) 1 ' 



i=l 



jS strin g(u) _ -q J R n+1 _ fejn+1 _ c . +1) (_2n + [i + A; - 1 - 2n] • 2 V ) . 



+ 

12... n 



(6.19) 

(6.20) 



k=l 



The reason why the conjugated P-matrices (jB8[) appear in this expression is that by commuting 
the <T 2 -matrices, arising from (j6.7j) . to the right, the relation 



R a b(u) = o z a R ab{u)<J z a = alR ah (u)a z b 



(6.21) 



is employed, cf. (|B7|) . 

SillCe [ P (l...n)' CT fn)] = and WtnY A ^...n) ^( n ) ^iL.n)] 

carries over to the fused ET~-matrices, 



0, the periodicity property (|6.5 



^ri2...n)( U + 7r ) 



-^<12...n»( U + 7r ) 



-!) n <n)^ 12 ... n )M<n) 

T) n C< n >-^<12...n>( n ) ,T <n> 



(6.22a) 
(6.22b) 



where the alternating sign results from successive application of (|2.5f[) . 



D. Fusion hierarchy for OBC 

From the fused quantities, it is again possible to derive a family of commuting operators 

rW(u) = st r<l ... n> { ^ <1 ...„»(n) T« 1 ... n »(n) P^ . n>> (u) f<i.„ ft> (« + [n - 1] • 2r?) } (6.23) 

that extends the existing familiy of commuting super transfer matrices t{u) = t^\u) such that 
[t^ (u) , (v)] = for all i,j > 1. The quantity P<i...n.>(it) appearing in (|6.23p is related to the 
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fused object 



r ( i... n) (« + [n - 1] • 2r?) = P+.„ ri(«)T 2 (ti + 2rj) ■ . . . ■ T n {u + [n - 1] • 2 V ) P+.„ 

= TT (6 24) 

C(«)C(« + 2»/) -...-CC" +[«-!]■ 2»/) • [ ■ ' 

in the usual way by restriction to the only relevant matrixblock after applying the respective A- 
transformation. In the case of general open boundaries the fusion hierarchy for n > 1 is found to 
be 

i \ / \ m/ . T (n+1) (u) A(u + \n - 11 • 2<n) , . , 1W N 

£ n (u) C(2w + 2n • 2r/) 

where A (it) labels the OBC super quantum determinant defined in (|D23P and 

n 

Zn(u) = []C(2n+ [n + k] -2?]) . (6.26) 

k=l 

The structure of this fusion hierarchy can be further simplified by introducing the rescaled quan- 
tities 

/n— 1 \ 



K{U) = C{2u + l 2ri) f(n){u) = " 1 II C 1 ^) ] t w (u) (6.27) 



A(«) 

,i=l / 

with the convenient definitions f (0) (-u) = — r (0) (u) = 1 and f (1) (u) = — r (1) (u) such that (|6.25p 
becomes 



f (n) (u) • f w (u + n-2rj) = f (n+1) {u) - A(u + [n - 1] ■ 2r/) • f (n_1) («) . (6.28) 

VII. TQ-EQUATIONS FOR OBC 

As in the PBC case, the fusion hierarchy (|6.28p provides a system of relations between the 

eigenvalues A^ n \u) of the fused (super) transfer matrices. Defining A(u) = A^(u) and after 
shifting u — > u — n ■ 2ij this yields 

A(n = r-r- ^ - A(u-2?? ^— v '1 7.1) 

A(«)(n-n-2r/) AW (u - n ■ 2rj) 

Introducing the functions 

h (n) ( n ) 7 (n) (n) g(n) (y) _ £(„) _ „ . 2r? ) (7.2) 
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where 

(n) ^ + 2,) » S1 n(2,-[ 2j -2].2,) 

' v ; sin(2n) sin ( 2n - [2j - 3] ■ 2rj) ' J 

n 

h^ n) (u) = -(-l) n Yl w+ sin ( n - k ■ 2rj - ■ ui~ sm(u - k ■ 2rj - (7.3b) 

fc=0 

the eigenvalues can be written as 

sin(2u) Q(" +1 )(n + 2r ? ) 



A(u) = R+(u)K s (u + 2r,) 



sm(2u + 2r)) Q( n )(u) 
A(u- 2r]) sm(2u - 2r]) qC"" 1 ) (u - 2??) 



(7.4) 



(« - 2rj)8T s (u) sin(2u - 4??) Q(n) ( u ) 

where the functions ^f(u) are defined in ()C7p , Now assume that there exists a limit Q(^) = 
lim ri _ >00 Q( n \u) which factorizes according to 

Q(u) = f N {u)q(u) with /(^^e 1 ^ 8 ^"^ . (7.5) 

Resubstituting A(u) = — A(u) by virtue of f)6.27[) it follows that 

a / \ ^ , . q{u-2r]) q{u + 2rf) 

A(u) = H a [u) — — ftn — — • ( 7 - 6 ) 

q(u) q[u) 



H«(u) = Sin g U + ^ £+(u - 2,)^(n) ' " ^ + ^ ' " 



where 



sin(2n + 2r/) a a \sin(ii + 2rj) sin(u — 2rj) 

T r. . sin(2u) > . . . . / — sin 2 (u) 

H (u) = 7-i — ^ &i(u)Si7(u + 2n) 7 — V r 

v ; sin(2u + 2r/) 5 w 5 v /y Vsin(u + 2rj) sm(u - 2r}) J 

such that H a (u)H s (u-2rj) = C~ 1 (2u)A(u-2jj). In the case q{u) = q(u) this is precisely the result 



(7.7) 



C9|) obtained by means of Bethe ansatz 



By matching the the asymptotically leading terms on the left and right hand side of their 
TQ-relations, the authors of [20] derived a set of constraints among the boundary parameters of 
the open (non-diagonal) XXZ Heisenberg spin chain which permitted them to employ Q-functions 
of the same factorized structure as in the case of diagonal boundary conditions, see Eq. ()4.9p . 
Unfortunately, this approach is not applicable to the small polaron model since in the non-diagonal 
case the leading term of H a,s (u) does not have the same power in z = e m as the leading term of 
A(u). However, it is possible to construct Q-functions which at least yield the correct behaviour 
in the asymptotic regime. Based on observations for small system sizes we suggest that in the case 
of non-diagonal boundary conditions the Q-functions can be written as 

q(u) = q(u) + p(u) ■ (fi+a. - a+/3_) (7.8) 
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with q(u) being the Q-function allowing for factorization (|08[) as in the diagonal case and another 
unknown function p(u) depending on 77 and the diagonal boundary parameters ip±. To determine 
p(u) the ansatz (|7.8p should be used in the TQ-equation (|7.6p together with the analytical properties 
of the transfer matrix eigenvalues, in particular their asymptotic behaviour (|6.14|) . 



VIII. TRUNCATION OF THE OBC FUSION HIERARCHY 

From here on, for the sake of readability, some of the functions introduced above will be equipped 
with a second parameter indicating for them to be taken at that particular value of the anisotropy 
j). For instance, K ± (u,p) = K ± (u)\ ri ^ p and so on and so forth. 



A. if-matrix truncation 



It is convenient to define the following functions 

sin(2r/ n ) 



±5{K ± (Tu-2 Vn ,r ]n )} 



2n 

n 



sin(2it + k ■ 2r] n 



TT^l-f) sin([n + l][u T ^])nn 



sin(2-u - 2 • 2rj n ) £j- sm(2rj n ) 

sin(2u +[i + j] ■ 2r) n ) 

pW) y-2) am ^ t ij L " + ™> ii ii — shns — • 

where = w ± (r/ n ) and to introduce the shorthand notations 

= <<n» • ^Ql...n»( U + 2r ?) 

The truncation identities for the boundary matrices can then be expressed as 
C<i... n > K^. x n:> (u,rj n -i) C^y n> 

-B«l...n-2»^ n _ 2 »( n ' 7 ?n-l)- B <l...n-2» * 



Mn-i(w) 



■1) 
•2) 



.3) 



■4) 
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B. OBC super transfer matrix truncation 

In order to be compatible with the truncation identities for the boundary matrices, the i?-matrix 
truncation identities (|5,2p need to be recast, this time employing the C transformation matrices, 



C<l...n> -R<l...n»g(w,7/n-l) C<i... 



V 



■Mn-liu) al 



C(u) of 74" 2 \u + 2r? n _i,77„_i) 



(8.5) 



where in slight contrast to definition (|5,3 



Kjp>(u,rj) = #<l...n> „n>g(«) £<<!... 



.6) 



such that for the single row monodromy matrix 



T^(u,rj) = C<i... n » ii<l...n> gjv (u) • . . . • i?<l...n>g 2 (' lt )- R <l-ri>?i(«) C<1 
= C<l...n> ?<l...n>(^) C<l... n > 



<l...n> 



(8.7) 



the truncation identity at rj = reads 

f[-K-i(«fnU<, 



(n-2). 



V 



Again it is convenient to introduce the C-transformed object 



f^(u, V ) = C<i...„» f« 1 ... n »(n) 



<l...n> 



J.10) 



to easily recognize the truncation identity 
f ln \u+[n-l]-2r ln - 1 ,ri n - 1 ] 



C N (u) C N (u + 2 % _i) . . . C*(« + [n - 1] • 2 Vn . 



.11) 



[WufniiKH 
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Now that the individual truncation identities for all the objects involved in the construction of the 
fused OBC super transfer matrix r^ n \u) are known, it can be shown by simple matrix multiplica- 
tion 5 of HO"), (EHD, (El - ) and ([HTTP that 

r (x+ \ 



T (n) (u^n-i) = str <1>->B> < 



x- 



X + - ^«i... n - 2 » { Y } + (-l) n X- (8.12) 



I V A 7 

with the placeholders X ± and y defined by 



X± = (±1) 



n-l 



[]r ff («+i!-2vi) 



.fc=0 



-M'-iH ^-i(«)^-iW ^-i(=F«)^i(±«) (8-13) 



and 



^= £«l...n-2» i^ 1 ... ri _ 2 >>( u + 2? M-l) CT<<n-2» fjl °»l X 

X T <1 ... n _ 2> (« + 2»/ n _ 1 ) CT<< n _2» ^<l...„-2»( n + 2r ?n-l) X 
X (n°«) ^<l-»-2>( U + 2 ^- X + K™ " 2 ) " ^ • 2r ?«-l) B <l...n- 



2> 



5.14) 



= 0n-l( U ) ^«l...n-2» ^<l... n -2>( U + 2, M-l) T «l...n-2»0 + 2 ?7n-l) X 

x K <i...n-2»( u + 2? 7n-i) f<i.„ n _ 2 >(u + 2r/ n „i + [(n - 2) - 1] • 2r7„_i) B^ 1 1 ... n _ 2> . 
In the second step of equation (|8.14|) relation (|D22() has been employed to get rid of the a z factors 
such that (|8.12p eventually yields the truncation identities for the OBC transfer matrices, 



(U, Th-i) =4>n-l(u) ■ 1 -<Pn-l(u) ' ^ 2) (« + 2 ^n-l , Vn-1 ) 



.15) 



where 0„ (u) and 4 > ni u ) are rather lengthy expressions given by 



\\C N {u + k-2r ln 



k=0 



Ml N (u) fi+(u)n n (u) [u+(-u)v n (u) + v+(u)u n (-u)] 



C(«) 



iV 



5.16) 



((u + n-2r) n 

In terms of the rescaled transfer matrices (|6.28p it is reasonable to introduce 



lid, 



u 



tier 1 !?) 

i=l 

iu ri w 



?7=»7n 



i=l 



</>n( u ) 



V=Vn 



n-2 



i=l 



5.17) 



V=Vn 



5 Due to the cyclic invariance of the supertrace, all the matrix objects in (|6.23[) may be conjugated by means of the 
C-transformation without changing the actual super transfer matrix. 
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which yield the respective rescaled truncation identities 

f( n V,Tfc_i) =&Li(u) ■ 1 •f ( "- 2) (n + 2r ?n _ 1 ,r ?n _ 1 ) . (8.18) 

IX. SUMMARY AND CONCLUSION 

Starting from structures provided by the Yang-Baxter algebra (|2.6p and the reflection algebra 
(|6.ip . (|6.2p we have set up the fusion hierarchies for the commuting transfer matrices r^ n \u) of the 
small polaron model with periodic and general open boundary conditions, respectively. Taking the 
limit n — > oo we have obtained TQ-equations for the eigenvalues of the transfer matrices which can 
be solved by functional Bethe ansatz methods for periodic and diagonal open boundary conditions. 
The resultin g sp ectrum coincides with what has been found previously using the algebraic Bethe 



ansatz 



p,y, 3 is 



For generic non-diagonal boundary conditions the U(l) symmetry of the model corresponding 
to particle number conservation is broken. Therefore, the algebraic approach cannot be applied as 
it uses the Fock vacuum as a reference state and relies on this being an eigenstate of the system. 
This situation is well known from the (ungraded) XXZ Heisenberg spin-1/2 chain with non-diagonal 
boundary fields where in spite of significant activities a practical solution of the eigenvalue problem 
for generic anisotropics and boundary fields is lacking. Here we have used the strategies employed 
for the XXZ chain to the (graded) small polaron chain: as for the spin chain the fusion hierarchy can 
be truncated at a finite order for anisotropies being roots of unity, rj = n/(2(p + 1)) [l^]. We have 
derived the corresponding truncation identities for all boundary conditions considered. For generic 
anisotropies the spectral problem needs to be studied within its formulation as a TQ-equation. 

To actually compute eigenvalues of the transfer matrices further steps have to be taken: for 
anisotropies being roots of unity the analysis of the truncated fusion hierarchy can be be done fol- 
lowing the steps which have been established for the XXZ chain 
on the boundary fields may arise. For generic anisotropies the situation is more complicated: in 
the ungraded XXZ chain a (factorized) Bethe ansatz for the Q-function given in terms of finitely 
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291 ] where additional constraints 



18-20 



2J. For 



many parameters was possible only if the boundary parameters satisfy a constraint 
graded models such a constraint may be absent: in the rational limit rj — > of the model consid- 
ered here the functional form of the Q-function remained unchanged when off-diagonal boundary 
fields where added 



261 ] . Similarly, the nilpotency of the off-diagonal boundary fields may allow 
for a general solution of the small polaron model. Starting with the proposed ansatz (|7.8p for the 
Q-function the derivation of Bethe type equations appears to be possible in the generic case. These 
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open questions shall be addressed in a future publication. 
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Appendix A: Graded vector spaces 

Fermionic lattice models exhibit a natural TL 2 gradation on their local space of states, i.e. 
V = Vo V\ is equipped with a notion of parity, 

p : Vi -> Z 2 , p( Vi ) ^ i G {0, 1} . (Al) 

Let dim Vo = m G N and dim V\ = n G N be finite. Then V is said to have dimension (m\n) and Vo, 
V\ are called the homogeneous subspaces of V . An element v G V is said to be even if p(v) = and 
is respectively called odd if p(y) = 1. While even elements of V correspond to bosonic states, odd 
elements represent fermionic states. For instance consider the case where both of the homogeneous 
subspaces Vb and V\ are one-dimensional such that the composite local space of states V = Vb © V\ 
is spanned by just one bosonic and one fermionic state. Then V is said to have -BF-grading, where 
BF refers to an ordered basis of V in which the first basis vector is associated to with the bosonic 
state (B) whereas the second basis vector is associated to the fermionic state (F). Now consider 
the tensor product of two copies of V. Taking into account the order of the basis states, the tensor 
product space will have BFFB-gr&d'mg, 

V (8) V = (Vo Vi) (8) (V © Vi) = (V © V ) © (V © Vi) © (Vi © V ) © (Vi © V x ) . (A2) 

B F F B 

In the following, the conventions from 



301 ] will essentially be adopted. 

Let {ei, e2, • • • , e m , e m+ i, . . . , e m+n } be a homogeneous basis of V, i.e. each basis element has 

distinct parity p(e a ), and for convenience let this basis be ordered, such that the first m elements 

span the even and the last n elements span the odd subspace of V, 

( 



p(a) = p(e a ) = < 



if 1 < a < m 

(A3) 

1 if m+l<a<m + n 



In order to deal with an algebra of linear operators, acting on the graded local space of states, it 
is necessary to extend the concept of parity to End(V), the space of endomorphisms of V. The 
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(m + n) x (m + n) basis elements of End(V) will be labeled e£ and are denned through their action 
on the above basis of V, 

efa = 5?e a . (A4) 
By extending the definition of the parity function to 

p(e/) = p(a) + p(/3) mod 2 (A5) 
End(V) becomes a Z2 graded vector space. A basis of the iV-fold product space 

End® N (V) = End(V) <g> End(V) ® . . . ® End(F) (A6) 

N v ' 

JV times 

can most naturally be obtained by embedding the local basis elements e a into this tensorproduct 
structure. Moreover, End® w (V) aquires a Z2 grading by a further extension of the definition of 
the parity function, 

P(ej? ® e/ 2 2 ® • • • ® ej 3 /) ee p(e + p(e^) + . . . + p(e//) mod 2 . (A7) 

When dealing with graded vector spaces, it is useful to replace the usual tensorproduct structure 
by a so-called super tensorproduct. The symbol ® s will be used to distinguish this new structure. 
With respect to a certain basis, the components of the super tensorproruct of two operators A £ 
End® k (V) and B £ End®*(F), where G N, is explicitly defined through 

(A ® s 5)% = (_1)IP(«)+P(fl] pM4« B\ . (A8) 

As pointed out in jscj, the super tensorproduct allows for a most convenient graded embedding of 
the e a into the j-th subspace of End® Ar (y), 

ej £ = t^-V 0s e P 0s ^(N-i) . (A9) 

A graded version of the permutation operator V is defined by the relation 

V(A ® s B) = (B ® s A)V . (A10) 

If (|A9|) is employed as a basis for End® (1^), the operator "P^ which permutes the i-th and the 
j-th subspace can explicitly be constructed as 

r ij = (-i)m e .j e .^ . ( aii) 

In the following, the definitions of some well-known operations, namely the matrix transposition 
and the trace operation, will be adapted to fit the needs of graded vector spaces. A nicely motivated 
and much more elaborated list of matrix operations on graded vector spaces can found in [161 ] . 



24 



Firstly, the super transposition of an element A S End(V) is denned by 

(A st ) a p = (_1)P(«)[P(«+/3)1^ . 



(A12) 



In contrast to the ungraded case, the super transposition is not an involution, i.e. applying 
the super transposition twice does not yield the identity operation. As pointed out in [8J, it 
is therefore convenient to introduce an inverse super transposition, 

(A ist ) a p = (-l)^^ a+ ^Ap a . (A13) 

The partial super transposition, i.e. a super transposition on the j-th subspace of End® N (V), 
is defined through 



(At ® s ...® s Aj® s ... ® s A N ) st > = Ai 



(A 



\st 



A N . 



(A14) 



The partial inverse super transposition is defined analogously. Please note that, as opposed to 
ordinary partial matrix transpositions on ungraded vector spaces, the successive application 
of partial super transpositions on all subspaces is gernerally not equal to a total super 
transposition, i.e. (A 1 <g) s A 2 ) stlSt2 / (Ai <g> s A 2 ) st • 



Secondly, the super trace of some A £ End(V) is given by 

str{A} = £(-!)*<«> A" 



(A15) 



For operators B G End® Ar (V r ) it is convenient to define a partial super trace on subspace j 
as 



ai ...Oj_i Qj+i ... a N 
%-l Pj+i ••• Av 



a N 

'3+1 ■■■0N 



(A16) 



Appendix B: Relation to Bracken's dual reflection algebra 



According to [8[ the dual reflection equation for quite general graded models reads 
R 12 (v - u)K+(u)R 21 (-u - v) ist ^ t2 K+(v) 

= K+(v)R 12 (-u - v) istlSt2 K+(u)R 21 (v - u) 



where 



i? 2 i(A) istlSt2 
^i 2 (A) istlSt2 



st 2 



isti 



(Bl) 



(B2) 
(B3) 
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By performing a super transposition on the first space and an inverse super transposition on the 
second, i.e. by applying (.) st i lst 2 t Q equation (jBf p one obtains the equivalent form 



R2i(v - u) stlist2 K+(u) s ^R 12 (-u - v)K+(vj lst2 

= K+(v) ist2 R 21 (-u - v)K+(u) stl R 12 (v - M ) stlist2 
In the case of the small polaron R-matrix as defined in (|2.3p one finds 

C(A) 



(B4) 



#21 (A) 



C(A-2t7) 
C(A) 



R 12 (X-4 V ) 
i?2i(A-4r?) 



(B5) 
(B6) 



C(A - 27?) 

At this point it is convenient to introduce a shorthand, which will henceforth be referred to as 
conjugated R-matrix, 



R ba (X) = M~ l R ba (X) M a 



(B7) 



with M being the so-called crossing matrix. For the small polaron model in particular, it is found 
that M = a z such that 



R ab {X) = R ba (A) = R ba 



(A) 



(B8) 



= Rl a W = <6 a (A) = R2W = CW • 

Using this conjugated R-matrix (jB8|) . the dual reflection equation may be written as 
R 12 {v - u)K+(u) stl R 21 (-u -v- Arj)K 2 + (v) ist2 

= K+(v) ist2 R 12 {-u -v- Ar])K+(u) stl R 21 (v - u) 

and is graphically depicted by 



^7, 



^ k;(u)«' 



(B9) 




K*(v)" 
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Appendix C: Algebraic Bethe ansatz for diagonal boundaries 



The reflection equation (|6.ip gives 16 fundamental commutation relations for the quantum space 
operators A, B, C and T> of which the following three are of particular interest, 

B(u)B(v) =B(v)B(u) (CI) 



v (v-u) s 2 (u+v) 
<Hv)s 2 (0) 

s 2 



; ) s 2(°) f s 9 (2v)s 2 (u+v) A(,,\ T>(iA\ 



(C2) 



T)(oAK(iA - s 4( u+v ) s 2( u - v ) K(~,\f)(nA s 2 (0)s 4 (2M)s (2t)) 



(C3) 



" s (u-v) s 2 (u+v) VKU/^KUJ (u) s 2 (2m) s 2 (u+v) s 2 (2v) 1 

using the abbreviation s fc (A) = sin(A + fcr/). To obtain the desired commutation relations, it is 
necessary to make an ansatz for a shifted P-operator 

D(A) = $(\)V(\) + <t>(\)A(\) (C4) 

and to determine the scalar functions 4>(X) and i?(A). It turns out that 0(A) = ;r^xy w hile $(A) 
remains arbitrary. Starting from the general boundary matrices given in (|6.3|) the diagonal case 
can easily be obtained by setting a± = f3 ± = 0. This leads to Bethe equations 

s 2^j) \ 2N _ s 2(^ ~ s 2(^ ~ VO TT s 4^j + s 2^j ~ ^) /^n 

s o(^) / SoC^i + ^+) s (uj + ^_) 11 s (z/j + v t ) s_ 2 (uj -ut) 

and super transfer matrix eigenvalues 6 

s 2(°) efrt..\ \ ( s a( n ) \ N TT s (u + ve) s 2 (ve - u) 



s 2(°) ^^ / s o( n ) A^-j^r s 4 (u + z^)s 2 (u-z^) 



51 ; V 5 s 2 (2n) Ql V Vs 2 (n)s 2 (-u); 11 



(C6) 



v s 2 (u)s 2 (-u)y Al s (u- u e )s 2 (u + v e ) ' 
Here ^^(m) label the diagonal entries (i.e. the eigenvalues) of the diagonal boundary matrices, 

fi~(u) = uj~ sin(^_ + u) , &+(u) = oj + sin(u + 2rj + ■0+) , 

(C7) 

&7 (u) = <jj~ sin(V>_ — u) , ^J(w) = w + sin(u + 2rj — ip + ) . 



Note, that this result corresponds to the one obtained by [13| . However, the authors of seem to have made a 
slight mistake when substituting their formula (57) into (61) to obtain (62), which should correctly read 

= sm(2tt + 4q) S in( U + t+) _ sin( M + 2 V - t+) ^ 
K 1 sin(2M + 2j7) v ; sin(2M + 277) v ' 
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Introducing the functions 



M 



l( u ) = Y\. sni ( n + 2?? + v£) sin(u — V() 



the eigenvalues f|CJ6[) can be recast as 



s 2 (2u) 



s 2 («) s 2 (-it) 



s 2 (2u)' 



S 3 («) s 2(~ n ) 



A? 



t/(u + 2ry) 



(C8) 



(C9) 



Appendix D: Super quantum determinants 



Consider a generic BFFB graded i?-matrix of the shape 

( a{u + 277) 





R{u) 



\ 



(Dl) 





a(«) a(2fj) 
a(2r?) a(u) 
\ -a(u + 2r])J 

where a(—u) = a(u) and o(0) = 0. At u = —2rj such an i?-matrix gives rise to a projector P~ onto 
a one-dimensional subspace 



2a(2?]) 



R(-2 V ) 



(° 





0^ 





1/2 


-1/2 





-1/2 


1/2 


v° 





0y 



Let T(u) be a representation of the graded YBA 

R 12 (u - v) Ti(u) T 2 {v) = T 2 (v) Ti(u) Rvziu - v) 
with the usual embeddings T\{u) = T(u) ® s 1 and T 2 (v) = 1 <g> s T(v), where 

V\(u) T\(u)\ ^ U(u) B(u) 
T\{u) T\{u)) bf ~ \C{u) D{u) i 

The PBC super quantum determinant (SQD) is defined as 



T{u) 



BF 



5{T(u)} = S ti 12 {P{ 2 T 1 (u)T 2 (u + 2rj)} 

= -{C(u)B(u + 2rj) - A(u)D(u + 2rj) 
-B(u)C(u + 2rj) - D(u)A(u + 2??)} . 



(D2) 



(D3) 



(D4) 



(D5) 
(D6) 
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At v = u + 2r\ and after dividing by a{2rf) the graded YBA yields the commutation relations 



C{u)B{u + 2r?) - A(u)D(u + 2rj) = C{u + 2rf)B(u) - D(u + 2r,)A(u) 
D(u)A(u + 2rj) + B{u)C{u + 2rj) = D(u + 2rj)A(u) - C(u + 2r])B(u) 
B(u)C(u + 2r?) + D(u)A(u + 2rj) = B(u + 2r/)C(«) + A(u + 2rj)D(u) 

These relations can be used to simplify the super quantum determinant to 

5{T{u)} = -[A{u)D{u + 2r?) - C{u)B{u + 2??)] . 



(D7) 
(D8) 
(D9) 



(D10) 



It remains to show that the super quantum determinant is a central element of the graded YBA, 
i.e. that it supercommutes with all the other elements A(v), B(v), C(v) and D(v) for arbitrary v. 
Consider the expression 



Ru(u - v)R 13 (u - w)R 2 3(v - w)Ti(u)T 2 (v)T 3 (w) . 



(Dll) 



Employing the graded YBE once it is obvious that 



v — > u + 2r\ 



(inil]) = R 2 3(v-w)R 13 (u-w) [22i2(u-u)2i(u)T 2 (u)] T 3 (w) 

■ -2a(2r])R2 3 (u -w + 2r])R 13 (u - w)P 12 T 1 (u)T 2 (u + 2r])T 3 (w) 



(D12) 
(D13) 



On the other hand, by applying the graded YBA relation twice it is found that 



v — > u + 2r\ 



doll]) = T 3 (w) [fli 2 (u-«)Ti(u)T 2 (t;)] R 13 (u - w)R 23 (v - w) 

■ -2a(2r ] )T 3 (w)P{ 2 T 1 (u)T 2 (u + 2rj)R 13 (u - w)R 23 (u -w + 2rj) 



(D14) 
(D15) 



Equating (|D13p and ()D15p and multiplying from both sides with P l2 gives 



(D16) 



{P{ 2 R 23 (u -w + 2 V )R 13 (u - w)P{ 2 }{P{ 2 T 1 (u)T 2 (u + 2 V )P{ 2 }T 3 (w) 

= r 3 H{P 1 - 2 ri(«)r a (ti + 2t,)P{ 2 }{P{ 2 r 13 {u - w)r 23 {u - w + 2 ?7 )p 12 } 

where additional P^ 2 projectors have been inserted by virtue of the appropriate triangularity con- 
ditions. After a change of basis to the P^ 2 eigenbasis via A± 2 as defined in (|3.14p . it is easy to check 
that application of the supertrace str 12 { . } yields 



<4 6{T(u)} T 3 (w) 
ai8{T(u)},T 3 (w) 



T 3 (w) a z 3 5{T(u)} 


. 



(D17) 
(D18) 
(D19) 
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Similarly one may introduce the object 

5 [f («)} = str 12 { P u f 2 (u) f x (u + 2 V ) } (D20) 



which obeys the exact same super commutation relations and by (|6.9p turns out to be proportional 
to the inverse of the above SQD. In particular for the considered iV-site small polaron model it is 
found that 

N 

5(u) = 5{T(u)} = -( N (u + 2 V ) ( D21a ) 

i=i 

1 N 

S(u) = 6{T(u)} = ( D21b ) 
* {U) i=i 

where qi labels the i-th quantum subspace (cf. section IV Bp . Moreover, the commutation relation 
(jD18[) extends to the fused quantities according to 

[a z <<n>> HTW},T <<L .. n>> H = (D22) 

with <r^ n> being defined in equation (|3.16p . 

In the open boundary case, the place of S {T(u)} is taken by another object A(u) which will most 
appropriately be called the OBC super quantum determinant. Generally, the SQD is what you get 
when you alter the first fusion step such that, instead of creating a higher dimensional transfer 
matrix by projection on a three dimensional auxiliary space, you now create a lower dimensional 
object by projecting onto the complementary one dimensional space. In a sense, loosely speaking, 
you do a reduction instead of a fusion and find that the open boundary SQD factors as follows, 

A(u) = str 12 { P 12 K+{u + 2rf)R X2 {-2u - %rj)K+{u)T x {u)R 12 {2u + 2 V )T 2 {u + 2rj) } 

= 5{K + (u)} -5{T(u)} -6{K-(u)} ■ 5{f(u)} (D23) 

' a ^^y5{K+(u)}.5{K-(u)} 

where T~ (u) was defined in (|6.8|) and 

5{K+(u)} = str 12 { P{ 2 K+(u + 2 V ) R 12 (-2u - 3 • 2rj) K+(u) } 

(D24a) 

= g(-2u - 6r?) • det {K+(u)} 

S{K'{u)} = str 12 { P{ 2 K{{u) R 21 {2u + 2rj) K 2 (u + 2rj) } 

'" ' " (D24b) 

= g{2u + 2rf) ■ det {K~(u + 2rj)} 

with the function g(u) being introduced in the context of (|2.5dp . Since a ± ■ f3 ± = 0, as mentioned 
in section [VI A\ the expressions det {K ± (u)} are well defined. 
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Appendix E: Transformation matrices 



This appendix presents a collection of matrix representations of the various similarity transfor- 
mations employed in this paper. It is convenient to define the coefficients 

-1/2 



a n = 



2n f\„] I r 1/2 and ft- f ^\v=m \ 
n + 1 WW™* and b -{[3] q \ v=m ) 



where [n] q denotes the usual g-deformation of an integer n G N defined by 

n _ -n 

[n]n = - V with q = ^ v 

q — q 1 

and to set 



A {1) = B 



<i> = C<i> 



1 
1 



(Ela) 



(Elb) 



(Elc) 



^(12) 

^<12> 
C<12> 



10 







1 1 







V2 V2 

1 







1 1 



o 



V " V2 V2 

diag(a 2 , l,a 2 ) 
diag(a 2 , 1, 1) 



(E2a) 

(E2b) 
(E2c) 



(123) 



1 








v° 

5«i2» = diag 
C«i23» = diag(a 3 , 1,1,1 















1 

V3 


1 

V3 





1 








1 

y/3 

















2V2 

3 


3 





V2 
3 


V2 
3 


2^ 
3 





V2 
3 








2^2 

3 











V2 
3 





03, 1, 1,03) 

























1 


1 











1 




















V2 
3 


V2 
3 





2y/2 

3 


V2 
3 






(E3a) 



(E3b) 
(E3c) 
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^4(1234) = — 





/ /ft 
' V2 



















1 

V2 


1 

V2 


(J 


1 

VI 




U 


u 


u 


1 

V3 


n 
U 









































3 
2 


1 

~2 












x 

— 2 


3 
2 







1 











5 
3 





77? 












3 
2 













I 

~ 3 
















1 

— 3 







u 


n 
U 


U 


u 


U 













1 

3 


o 













1 

3 
























I o 














diag(a 4 , 1,6, l,a 4 ) 

















-E><1234> 

C«i234» = diag(a 4 ,l,6, 1,1) 



l l 

V3 y/3 



























_ 1 

3 













1 

V2 









3 
2 








1 

V2 



o \ 







_! 

2 

_! 

2 



_ 1 

2 
































1 

V3 


1 

V3 





1 

V3 

















1 

V2 





1 

V2 


1 

V2 























V2 












































1 

3 


1 

3 





1 

3 
































1 

3 


1 

3 





1 

3 











1 

3 


1 

3 





1 

3 





















5 
3 









3 
2 



"3- 

4 o 








_ 1 

2 




_ 1 

2 

_ 1 

2 








o / 

(E4b) 

(E4c) 
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